Theory of Hysteresis Loop in Ferromagnets 
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A theory of the hysteresis loop in ferromagnets controlled by the domain wall motion is presented. 
Domain walls are considered as plane or linear interfaces moving in a random medium under the 
action of the external ac magnetic field H = Hq smujt. We introduce important characteristics of the 
hysteresis loop, such as dynamic threshold fields, reversal field etc. together with well known char- 
acteristics as coercive field and hysteresis loop area (HLA) A. We show that all these characteristics 
are regulated by two dimensionless combinations of the Ho and lu and intrinsic characteristics of the 
ferromagnet. The moving domain wall can create magnetic bubbles playing the role of pre-existing 
nuclei of the reversed magnetization. We discuss a simple model of this process. For magnetization 
reversal determined by domain inflation we predict that HLA scales as ^ oc uj^Hq with a = 1/2 
and /3 = 1/2. Numerical simulations confirm this result. 
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Magnetic hysteresis loops (HL) have been first studied 
already more than a century ago 0]. However, the un- 
derstanding of this process in thin magnetic films is still 
rather poor. Many efforts have been devoted recently to 
the prediction [|[j5) and experimental verification of 
the scaling behavior of the hysteresis loop area (HLA) as 
a function of the frequency and amplitude of the applied 
magnetic field in thin magnetic films The scaling 
behavior of the HLA has been first reported in the pio- 
neering work § for 3D magnets. While there exists an 
extended literature on the hysteresis of 3D magnets, the 
properties of the HL in 2D systems are much less known 
P-p^ . Critical exponents found in the experiments with 
thin films differ dramatically for different materials |^^] 
and probably for different regimes. Different authors dis- 
agree with each other |^|-^ and also disagree with numer- 
ical simulations j^]. 

More recently mean field type models with single ||] 
or many pO| ] relaxation times have been applied to ana- 
lyze the experimental data for the description of the HL 
controlled by nucleation processes. These authors pre- 
dict a logarithmic dependence of the coercive field He 
on the rate of the applied magnetic field H. In a, recent 
experiment it was found that the HLA depends on 
frequency lo of the applied field as power with a small 
exponent a (~ 0.03 — 0.06) or, possibly, logarithmically. 
In the framework of the same approach the HLA must 
behave also logarithmically in H. However, such a de- 
pendence has never been observed experimentally. 

In this article we propose a new analysis of the HL 
by formulating a rather general approach to the magne- 
tization reversal mechanisms. In particular, we indicate 
several important measurable characteristics of the HL 
besides of the HLA. It turns out, that these characteris- 
tics are governed by two dimensionless parameters com- 
bined from the field frequency lo, its amplitude Hq and 



characteristics of the magnetic material. Everywhere in 
what follows we assume that the external field varies har- 
monically in time, H{t) — HosiriLut. 

In general, hysteresis behavior may have various ori- 
gins. It can be mediated by nucleation processes, by the 
domain wall (DW) propagation or simply by retardation 
of the magnetization due to fluctuations. In this article 
we will restrict ourselves to the second mechanism, i.e. 
the hysteretic behavior due to the finite velocity of the 
DW propagation, and establish conditions at which it is 
dominant. Discussion of other mechanisms as well as a 
more detailed presentation of the results of this paper 
will be postponed to a future publication 

We consider magnets of Ising (uniaxial) symmetry. 
Their properties may be very different depending of the 
strength of the anisotropy. In the experimentally stud- 
ied films the anisotropy was very weak. In this case the 
domain wall width is large in comparison to the lattice 
constant. On the contrary, in the original Ising model 
the anisotropy is assumed to be large and DW width I 
is simply the lattice constant. However, these different 
models becomes equivalent after a simple rescaling; the 
DW width should be accepted as a new elementary (cut- 
off) length. It means that we consider a spin cluster of 
the linear size I as a new elementary spin. 

Disorder plays an important role in the DW propa- 
gation: it creates a finite threshold value Hp the exter- 
nal field has to overcome to move the DW and lowers 
the DW velocity considerably in the low-field regime, 
{H — Hp) <C Hp, where the DW motion shows critical 
behavior ||l2|,|l3| . Peculiarities of the two-dimensional sit- 
uation are: much higher mobility of the DW as well as 
much stronger fluctuations. This makes the experimental 
situations much more diverse than those for a 3D magnet. 

We first consider the individual DW motion by for- 
mulating an effective equation of motion under the in- 
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fluence of the field H{t), which we solve then in a finite 
geometry. This will lead us to the characteristic fields 
Hti, Ht2, He, Hr and the HLA A which we analyze in 
several limiting regimes in which simple power scaling is 
valid. 

In general, the DW motion in an impure magnet is 
highly non- linear. However, as it was shown in | p2| , p^ , 
after integrating out DW fluctuations on time scales less 
than the dynamical correlation time t^ defined below, 
the effective equation of motion for the center of mass 
coordinate of the interface Z is given by 



/ H 



X , X ^ 1 



(1) 



where 7 represents the bare DW mobility. Inside the 
critical region H — Hp ^ Hp, f{x) obeys a power law 
with an exponent of the order 1/3 from an e-expansion 
around five dimensions extrapolated to to d~2 dimen- 
sions ||l2j. On time scales t ^ ty (i.e. for uity <C 1), 
and length scales L ^ fy where ty ~ —jriir — 1) 

lip lip 

and ^„ oc [jj- — l)^" , fluctuations of the DW around its 
mean position are weak and will be neglected therefore 
in the following. Here u and u are the correlation length 
and time, respectively, exponent of the driven interface, 
with v{d = 2) = 1 and v{d ~ 2) ~ 4/3 again from an 
e-expansion around five dimensions. Outside the critical 
region, H 3> Hp, 9 = \ and the domain wall is macro- 
scopically flat. 

We solve now equation of motion for a planar DW 
which can move over a distance L forth and back. Here 
L corresponds to a characteristic size of the system. In 
the simplest case L denotes its linear extension. In some 
experiments only one domain wall survives and the model 
problem of rectilinear domain wall motion is close to re- 
ality. Another situation may occur in systems with rare, 
but strong extended defects. In passing such a defect do- 
main walls may form bubbles of reversed spins. These 
bubbles play then the role of prepared nuclei in the next 
half-cycle of the magnetization reversal. In this case L 
can be identified with the mean distance Led between the 
extended defects. E.g. in random field magnets strong 
defects result from a coherent fiuctuation of random fields 
on neighboring sites in which case L^d is given by 



L 
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(2) 



where rj > Hq and T denote the strength of the uncor- 
related random field and the DW stiffness, respectively. 
Here T/l ^ r/ ^ Hp and c is a constant of order unity. 

The domain wall is assumed to be fixed at the left 
boundary of the sample Z = at the initial moment. 
Solving equation of motion (|l|) for the domain wall co- 
ordinate Z in harmonically oscillating magnetic field, we 
replace the integration over time by the integration over 



the field, using dt ^ ^ 



dH 



This yields 



Z{H) = 



iHp 



^ Jh. 



H 



H — Hr, 



dH 



Hp J ^Hl - m 



(3) 



This equation is correct for H > Hp. At smaller val- 
ues of H the domain wall does not move. The sign of 
the square root changes each time as H reaches its max- 
imum or minimum value ±iJo- The second necessary 
prescription is to substitute H — Hp by —{H — Hp) in the 
argument of the function / and Hp by —Hp in the lower 
limit of the integral (||) when H is negative. To transfer 
from the coordinate Z to the magnetic moment AA, we 
write M.{Z) = Msi^Z — L)/L. The area of hysteresis 
loop A can now be expressed in integral form: 



A = AHpMs{l 



1 



LH 



Z{H) dH ). 



(4) 



p Jh- 



For later applications, it is useful to introduce the func- 
tions $(a, 6) and ^(a) defined by 



Ha.h)^ fix -I) 



dx 



(5) 



and $(a) = $(a, a) . Eq. (||) can then be rewritten in 
the form Z/L — u^^^{H / Hp,v), where u and v are the 
two dimcnsionlcss variables 



= LoL/{-lHp), v = Ho/Hp 



(6) 



It is clear, that the HLA A also depends on these two 
parameters. 

There are two important characteristic values for the 
amplitude Hq which separate hysteresis loops of differ- 
ent shapes. The first of them is the dynamic threshold 
field Hti, which is the smallest value of Hq at which the 
domain wall reaches the right boundary of the sample 
Z = L. At Hq < Hti the magnetization is not reversed 
fully, the hysteresis loop is asymmetric, whereas at larger 
values of Hq the hysteresis loop is symmetric under in- 
version H — > —H,M — > —M (Fig. 1). The value of 
Hti is determined by solution of equation 



2-fHp 



H„ 



(7) 



At -ffo = Hti the hysteresis loop is symmetric with re- 
spect to reffections in the axis H and A4. The second 
threshold field Ht2 is defined as a value of Hq at which the 
domain wall reaches the right end of the sample Z = L 
during one fourth of period, just aX H = Hq. The def- 
inition of Ht2 differs from Eq. (^) by the absence of 
factor 2 in denominator of the l.h.s. The hysteresis loops 
corresponding to Hq > Ht2 acquire characteristic "mus- 
taches", single- valued pieces of the curve Ai{H), which 
are absent in hysteresis curves for Hq < Ht2 (see Fig. 
Id). 



2 



At a fixed Hq > Hti it is possible to define the coercive 
field He by the requirement A4{Hc) = 0. The so-called 
reversal field is defined as the field value at which the 
magnetic moment reverses fully. At H = Hr the two 
branches of the hysteresis curve intersect each other, i.e. 
Z{Hr) = L. At H between Hr and Hq the magnetic mo- 
ment remains a constant A4 = A4s- The values He and 
Hr are shown in Fig. lc,d. Using equation (||) and the 
relation for M{z), we find for H^- 



2jHp 



*(—,—)• 
Hp Hp 



(8) 



For Hf2 < Hq the relation for Hr differs from (||) by 
the absence of the the factor 1/2 in the l.h.s. For 
Hti < Hq < Ht2 the reversal field Hr is given by 



' H„ 



Hp Hp 



(9) 



Thus, also the ratios H^/Hp and H^/Hp are functions of 
u and V only. 

We analyze different asymptotics of the hysteresis loop 
characteristics starting with the simplest limiting case 
u ^ 1. In this case the contribution to Z{H) from 
the critical region is negligible and no further restric- 
tion for the applicability of (|l|) apply. The main con- 
tribution to the integral (||) comes from large values 
of the argument and the function $(a, 6) can be calcu- 
lated explicitly by replacing f{x — 1) by x and setting 
the lower cut-off for the integration equal to zero. This 
yields $(a, 6) w 6 



As a result Hp disappears 



from all expressions and a new dimensionless parameter 
u — LoL/ {^Hq) emerges. In particular we get for dynamic 
threshold field Hti'> Hp 



Hti 



27 



(10) 



Similarly, He and Hr are determined by He = 
^Hn{2Ho~Hti) and Hr = ^ Ht2{2Ho - Ht2). The 
shape of the HL for Hq > Ht2 is described by 



M = Ms 1 - 2u(l - - {H/HoY) 



(11) 



It can be used also for the range of amplitudes Hti < 
Hq < Ht2 on the lower branch of the hysteresis curve 
< H < Hq. On the upper branch of this curve 
Hr < H < Hq the sign of the square root must be re- 
versed. If Hq > Ht2, the HLA is determined as follows: 



A 



AMsHr 



1-u 



Ho 



(12) 



On the contrary, for Hq < Ht2, the area is: 



A = 4MsHr (1 



TtMsHqU. 



(13) 



In the range of existence of the full HL u < 1/2. There- 
fore u cannot be large. However it can be very small, i.e. 
ujL ^ ^Hq. In this case we find Hr ~ He = Hq^/U = 
A/wiiJo/7, M{H) ^ Ms{l - 2H^/{uH^)) and 



A W iMsHr OC W^/^^y^ 



(14) 



Thus, for -u ^ 1 we find the scaling behavior of all hys- 
teresis characteristics with universal critical exponents, 
independent on the pinning centers. Therefore, one can 
expect that the same scaling is valid for a clean ferro- 
magnet with the relaxational dynamics. 

Next we study the opposite limit of small w <C 1. In 
this case Hti and Ht2 are close to Hp as follows from 
(0)' (^' ®- Solving the integral (^) approximately and 
employing for f{x) at small x (^, we get 



f&-l)''+i/2 1 a-1 



(15) 



where we introduced the incomplete _B-function defined 
by the integral 



B{x,y,z) = [ s^'-^il- sf'^ds 
Jo 



(16) 



In particular, we get $(a) « (a - l)'^+^/ %[4+t/2) ■ 
One then finds from (^ for Hti 



Hti 

Hn 



(1 + l/26')r(0 + 1/2) ujL 



27rr(6l) iHp 



(17) 



Ht2 is given by the same expression if we replace L by 
2L in (|l^. The ratio {Ht2 - Hp) /{Hti - Hp) = 22/(2^+1) 
does not depend on u. We omit here the explicit expres- 
sions for He and Hr but give finally the results for the 
HLA in the case u <C 1. For (i) — 1 <C 1 we get: 



A 



iMsHp 



1-f^— ^ '- / B{e + l,l/2;x)dx. 



(18) 

where w = {Hr — Hp)/{HQ — Hp). In the case (ii) 
w ^ 1, uu <C 1 we find: 



A w 4MsHp 



{{0 + l)uv)'/^'+'^ 



(19) 



Finally, for (iii) uv ^ 1, Hr ^ Hp and therefore all es- 
sential results formally coincide with those for u ^ I . 
In the cases (i) and (ii) A is close to = 4:A4 sHp but 
the deviation A— Aq scales with the parameters u and v. 
The prerequisite inequalities uty <C 1 and <C L, lead 
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for u ^ 1 to further restrictions on L and w, which in 
general can be fulfilled if L is large enough . 

We have also used Monte-Carlo simulation with 
Glauber dynamics to check (|l|),(|l^) and ( p^ for a ran- 
dom bond Ising model. The disorder has weak influence 
on phase diagram, but decreases significantly the DW 
velocity. However, no measurable pinning threshold field 
Hp has been detected. To separate domain growth from 
domain nucleation we have studied two different cases: 
(i) specially prepared defects for fast nucleation, (ii) gen- 
erated nucleus with the opposite direction of magnetiza- 
tion in the form of circle or stripes. The magnetization 
changed according to the model of a straight DW, i.e. 
similar to Eq. (|ll|). In this case usually one or two do- 
mains inflate from the nucleation centers. The distance 
L jv between nucleation centers play the role similar to the 
system size L in the above consideration. We have stud- 
ied scaling behavior of HLA A (see Eq. |lj) for the case 
when it is determined by the domain growth. Details will 
be pubfished elsewhere [0. We found scaling A (x cu^ 
with p = 0.49 ± 0.05 over three decades of co. However 
the range of Hq values available is only one decade. For 
this reason we have only checked that the dependence is 

1 /2 

consistent with A oc Hq . 

To conclude, we have studied the hysteresis process in 
controlled by the DW motion. We introduced two dy- 
namical threshold fields Hti and Ht2 corresponding to 
the occurrence of the full magnetization reversal and to 
the occurrence of the single-valued parts on the hystere- 
sis curve respectively. These dynamical threshold fields 
are larger than the static threshold field Hp which is re- 
quired to start the motion of the DW. We established 
that Hti and Ht2 expressed in units Hp are functions of 
one dimensionless parameter u = wL/^Hp . The coer- 
cive field He and the reversal field Hr expressed in the 
same units are function of two dimensionless parameters 
He - HpF (5^, f ) and H^ = HpF (^, f ) . Ex- 
perimental observation of this type of scaling would be 
the best indirect evidence of the DW motion controlled 
hysteresis. A direct observation of the DW motion in 
principle is possible [Q. At large fields Hq ^ Hp the 
defects are inessential. Therefore the dependence on 
Hp must vanish from all scaling laws. It happens in- 
deed, and both fields He and Hr are expressed in terms 
of one dimensionless parameter loL/^/Hq We presented 
also corresponding equations for the HLA to which the 
most experimental efforts were concentrated. However, 
we would like to emphasize that the HLA is not the 
only measurable characteristics of the HL and even not 
the most informative of its characteristics: the fields 
Hfi, Ht2, He, Hr as well as the shape of the hysteresis 
curve are not less interesting. In the case when the driven 
DW are almost free {Hq S> Hp) and the HL is narrow 
{Hq ^ Hti) the HLA was found to be proportional to 
lo^^'^Hq^'^. This conclusion is supported by our numerical 



MC simulation. 
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FIG. 1. Schematic pictures of hysteresis loops. 

(a) Incomplete HL for Hq < Hti. 

(b) Symmetric HL for Ho = Hti. 

(c) The HL for Hti < H„ < Ht2- 

(d) The HL for Ho > Hr2. 

The values Hp, He, Hr, Ho are marked in all figures. 
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